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ABSTRACT 

We examine the dynamics of stellar systems embedded within cold dark matter (CDM) halos in 
order to assess observational constraints on the dark matter content of Local Group dwarf spheroidals 
(dSphs). Approximating the stellar and dark components by King and NFW models, respectively, 
we identify the parameters of dark halos consistent with the kinematics and spatial distribution of 
stars in dSphs as well as with cosmological N-body simulations. Our analysis shows, in agreement 
with previous work, that the total mass within the luminous radius is reasonably well constrained 
and approximately independent of the luminosity of the dwarf, highlighting the poor correspondence 
between luminosity and halo mass at the extremely faint end of the luminosity function. This result 
implies that the average density of dark matter is substantially higher in physically small systems 
such as Draco and Sculptor than in larger systems such as Fornax. Because massive CDM halos 
are denser than low mass ones at all radii, these results imply that Draco formed in a halo 5 times 
more massive than Fornax's despite being roughly 70 times fainter. Stellar velocity dispersion profiles 
(<7 p (i?)) provide further constraints; in systems where data exist, <r p (R) remains flat almost to the 
nominal "tidal" radius, implying that stars are deeply embedded within their cold dark matter halos 
and are therefore quite resilient to tidal disruption. We estimate that halos would need to lose more 
than 90% of their original mass before tides begin affecting the kinematics of stars, but even then the 
peak circular velocity of the dark halo, V max , would be barely affected. We estimate that V max is about 
3 times higher than the central velocity dispersion of the stars, a result in agreement with previous 
estimates and that alleviates significantly the CDM "substructure crisis". We use these results to 
interpret the structural differences between the M31 and Milky Way (MW) dSph population and, in 
particular, the observation that M3J dwarfs are physically more extended by approximately a factor 
two than MW dwarfs of similar luminosity. Our modeling indicates that the difference in size should 
be reflected in their kinematics, and predicts that M3J dwarfs should have velocity dispersions up 
to a factor of ~ 2 higher than their MW counterparts. This is an eminently falsifiable prediction of 
CDM-motivated models of dSphs that may be verified with present observational capabilities. 
Subject headings: 



1. INTRODUCTION 

Ever since the first measurements of their velocity dis- 
persion became available (Aaronson J983, Aaronson & 
Olszewski 1987), dwarf galaxies in the vicinity of the 
Milky Way have been a difficult puzzle to piece together 
in models of galaxy formation. Their relatively large size, 
low luminosity, and sizable velocity dispersion suggest 
the presence of large amounts of dark matter (Arman- 
droff, Olszewski & Pryor 1995, see the review of Mateo 
J998; for more recent work consult Kleyna et al. 2005, 
Munoz et al. 2005, 2006 and references therein). How- 
ever, the absence of clear correlations between inferred 
dark matter content and the structural properties of the 
luminous component have led to arguments about the 
true physical nature of these objects and to various pro- 
posals to explain their large mass-to- light ratios. 

Velocity dispersions unduly affected by binary stars 
(Olszewski et al f996); line-of-sight alignment of un- 
bound stars (Kroupa J997), and remnants of stellar clus- 
ters severely disrupted by tides (Metz & Kroupa 2007) 
have all been considered in the literature, but consensus 
now seems to have been reached. Dwarf spheroidal galax- 
ies (dSphs, for short) are widely regarded as dark matter- 
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dominated systems sampling the extreme low-mass end 
of the dark halo mass function. Nevertheless, a number 
of issues regarding the dark matter content of these elu- 
sive systems, as well as the spatial extent of their dark 
halos, remain largely unresolved. 

These issues were brought into focus when cosmolog- 
ical N-body simulations revealed the presence of sub- 
stantial substructure in galaxy-sized cold dark matter 
(CDM) halos (Klypin et al. 1999, Moore et al. 1999). 
These simulations indicate that hundreds of self-bound 
cold dark matter halos massive enough (in principle) to 
harbor a dwarf galaxy are expected to populate the halo 
of the Milky Way, in sharp contrast with the mere tens 
of dwarf galaxies known to orbit the Galaxy. This real- 
ization rekindled interest in explaining the detailed cor- 
respondence between dark halos and luminous galaxies 
at the faint end of the luminosity function, an issue that 
had long been highlighted as a challenge for hierarchical 
galaxy formation models (White & Rees 1978, Kauff- 
mann, White & Guiderdoni 1993). 

The leading scenario for reconciling the discrepancy be- 
tween "luminous" and "dark" substructure in the Milky 
Way halo envisions dwarf galaxies as able to form only in 
halos above a certain mass threshold. The threshold is 
determined by the need to retain gas and to sustain con- 
tinuing star formation despite the effects of feedback from 
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evolving stars and the heating from photoionizing radi- 
ation (Efstathiou 1992, Bullock, Kravtsov & Weinberg 
2000, Somerville 2002, Benson et al. 2002). Adjusting 
the mass threshold appropriately, the scarcity of lumi- 
nous dwarfs may be explained by the relatively few mas- 
sive substructure halos that exceed the threshold (Stoehr 
et al 2002, Hayashi et al 2003, Kazantzidis et al 2004). 
This is an appealing and elegant solution, and offers a 
relatively clean prediction that may be tested observa- 
tionally: most dwarfs at the extreme faint end of the 
luminosity function should inhabit relatively massive ha- 
los, of mass comparable to that defined by the threshold. 

For example, given the large spread in luminosity of 
dwarfs in the Local Group 3 , and the expectation that 
they should all inhabit halos of similar mass, this pro- 
posal implies that there should be little correlation be- 
tween the dark matter and luminous content of a dwarf. 
Massive dark halos would also be more affected by dy- 
namical friction (Penarrubia et al. 2002, Zentncr & Bul- 
lock 2003, Penarrubia & Benson 2005, Penarrubia et al. 
2006), leading to a bias in the spatial distribution of 
dSphs within the halo of the Milky Way relative to the 
bulk of substructure halos — an effect that may be poten- 
tially observable. 

Finally, their relatively massive halos would make these 
extremely faint galaxies more resilient to tidal disrup- 
tion, suggesting that tidal tails and other obvious evi- 
dence of tidal stirring should be present only in systems 
that have lost most of their original mass by the action 
of tidal forces. In these models, the "tidal" radii at- 
tached by King-model fitting to the surface brightness 
profiles of dSphs, as well as the "break" radii identified 
in the outskirts of some dSphs, should reflect just the 
edge of the luminous component rather than a feature 
of dynamical significance (e.g. UMi, Martmez-Delgado 
et al. 2001, Palma et al. 2003; Carina, Majewski et al. 
2005). Establishing the present dark matter content of 
dwarf galaxies in the Local Group thus promises to be a 
fruitful enterprise. 

Validating (or refuting) these theoretical expectations 
through observation is, however, not straightforward. 
One reason is that few dwarfs have gas on circular orbits 
and, therefore, stars are the only viable dynamical tracer. 
The interpretation of the observations is thus compli- 
cated by degeneracies between orbital shapes, their ra- 
dial dependence, and the overall mass profile (Wilkinson 
et al. 2002, Kazantzidis et al. 2004, Mashchenko et al. 
2006, Kleyna et 2004, Tolstoy et al. 2004, Wilkinson et 
al. 2004, Muhoz et al. 2005, 2006, Wang et al. 2005, 
Walker et al. 2006a,b, Sohn et al. 2006). 

The second reason concerns the fact that these dy- 
namical tracers lie, by definition, within the luminous 
radius of the dwarfs. Since dark halos are expected to 
be extended objects reaching far beyond the luminous 
confines of a galaxy, a certain amount of uncertain ex- 
trapolation appears inevitable. In this respect, mass- 
follows-light models (Richstone & Tremaine 1986) only 
provide an approximate estimate of the amount of dark 
matter within the luminous radius that cannot be safely 

3 We consider in this paper the "traditional" dwarfs brighter 
than M v ~ — 8 but recent discoveries based on the SDSS have un- 
covered the presence of many fainter dwarfs, extending the galaxy 
luminosity function to as faint as a few thousand Lq (e.g. Irwin 
et al. 2007 and references therein) 
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Observational properties of the Local Group Dwarf 
Spheroidals considered in this paper 
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extrapolated to larger distances. Finally, Galactic tides 
may affect dark matter and stars differently, especially 
if stars are, as expected, strongly segregated relative to 
the dark halo. Thus a dwarf may today inhabit a rela- 
tively low-mass halo (say, below the threshold mentioned 
above) even though it actually formed in a massive one 
that has since seen much of its mass stripped away by 
tides (Kravtsov, Gnedin & Klypin 2004). 

This paper is the first in a series that attempt to ad- 
dress these issues by interpreting available observational 
data on Local Group dwarfs in the cosmological context 
defined by the leading paradigm of structure formation; 
the cold dark matter theory. We start by considering 
the constraints on CDM halos placed by the structural 
and dynamical properties of the dwarfs, and use these 
results to try and interpret the origin of structural dif- 
ferences in the population of dwarfs that orbit M31 and 
the Milky Way, respectively. A future paper on the sub- 
ject will analyze in detail the effects of Galactic tides, 
both in the survival of the stellar components, as well as 
on the observational signatures imprinted by stripping 
(Penarrubia, Navarro & McConnachie in prep.). 

The paper is organized as follows. Section 2 introduces 
the relevant observations and the modeling procedure, 
whereas §3 applies the models to the Milky Way dwarfs. 
We apply these results to M31 dwarfs in §4, and conclude 
with a brief summary in §5. 

2. PRELIMINARIES 
2.1. Summary of Observations 

Table [T] lists the observational parameters of the Local 
Group dwarfs that we study here. We split the dwarfs 
into two groups, and list each in order of decreasing lu- 
minosity. The first set of objects are Milky Way (MW) 
dwarf spheroidals (dSphs), whereas the second set are all 
dSph type systems for which accurate core and tidal radii 
from King-model fitting are available. Most of the latter 
satellites orbit around M31, except for Cetus which is 
one of only two isolated dSphs in the Local Group. 

Table [1] lists the Galactic coordinates of each dwarf 
(I and 6), as well as its heliocentric distance, D. Dis- 
tances are taken from McConnachie et al (2004, 2005) 
where available, otherwise from the compilation by Ma- 
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teo (1998). Absolute magnitudes (M v ), core (R c ) and 
tidal (Rt) radii for the best-fitting King profiles to the 
surface brightness distribution are taken from Irwin & 
Hatzidimitriou (1995) for the MW satellites and from 
McConnachie & Irwin (2006) for Cetus and the M31 pop- 
ulation, and are quoted after rescaling to the adopted 
distance. One exception is And IX, with core and tidal 
radii taken from Harbeck et al. (2005), distance from 
McConnachie et al. (2005), and the absolute magnitude 
from Zucker et al. (2004). The other is And II, the 
one dSph in the Local Group whose surface brightness 
profile shows evidence for the presence of more than one 
dynamical component (McConnachie & Irwin 2006, Mc- 
Connachie, Arimoto & Irwin 2007) and where single King 
model fits are a poor description of the structure of the 
dwarf. In this case, we quote the And II core radius used 
by Cote et al. (1999), since this is consistent with the 
traditional definition of core radius as the distance from 
the center where the surface brightness drops by a factor 
of two. 

The central velocity dispersions (c p [0]) for the MW 
population in Table [1] are taken from the compilation 
by Mateo (1998) (his Table 7), although we note that 
since the publication of that review velocities for hun- 
dreds of stars in these dSphs have been measured with 
the aid of multi-object spectrographs. These data con- 
strain the velocity dispersion profiles, a p (R), in these 
systems, which are found to be approximately flat out 
to the nominal tidal radius (see, e.g., Fornax, Walker et 
al. 2006a, Battaglia et al. 2006; Leo I, Koch et al. 2006; 
Sculptor, Tolstoy et al. 2004, Westfall et al. 2006; Leo II, 
Sohn et al. 2006; Sextans, Walker et al. 2006b; Carina, 
Munoz et al. 2006; Draco and Ursa Minor, Wilkinson et 
al. 2004, Muhoz et al. 2005). Note that the presence 
of several distinct components may affect the measured 
kinematics of a dwarf, depending on the spatial distri- 
bution of the tracers for which velocities are available 
(see, e.g., McConnachie, Peharrubia& Navarro 2006, and 
references therein). We neglect here complications that 
arise from this issue, although we plan to address this in 
future work. 

Kinematic data for the M31 dwarf population is scarce 
and, comparatively, of poorer quality. For example, 
Chapman et al (2005) find that the velocity dispersion 
of And IX is ~ 6.8 ± 3 km/s, but note that adding a 
single (possible member) star to their sample raises this 
estimate to ~ 12 km/s. Similarly, the velocity dispersion 
for And II (Cote et al. 1999) is based on just seven stars 
and is thus subject to sizable uncertainty. The a p (0) es- 
timate for Cetus (Lewis et al 2006) is based on a larger 
sample of stars and it is thus more reliable. We present 
these data in Table [T] between parenthesis in order to 
emphasize that these data are not of the same quality as 
is available for the MW dSphs. We return to this issue 
in §4. 

The past few years have seen a dramatic increase in 
the number of dwarf satellites discovered around M31 
and the MW, in particular as a result of the completion 
of the Sloan Digital Sky Survey (Zucker et al. 2004, 
2006a,b,c, Willman et al. 2005, Belokurov et al. 2006, 
2007, Martin et al. 2006, Irwin et al. 2007). Many of 
these systems are morphologically not unlike the ones 
listed in Table [1] but they are typically much fainter and 
of much lower surface brightness, which has precluded 



robust fits to their surface brightness profiles. 

The sample of Milky Way dwarf galaxies used here only 
includes (i) dwarf spheroidal galaxies (i.e., non-rotating 
systems with little or no detectable gas). These galaxies 
must have: (i) King profile fit parameters available from 
the literature and (ii) measured central velocity disper- 
sion. Thus our sample excludes many of the recently dis- 
covered dSph galaxies in the Milky Way and M31. We 
have also excluded the Canis Major dSph (Martin et al. 
2004) because of its uncertain nature. A lively debate 
can be found in the literature, with arguments in favour 
of this system being a new galaxy (Martmez-Delgado et 
al. 2005, Bellazzini et al. 2006) and against it (e.g. Mo- 
many et al. 2006, Moitinho et al. 2006). Sagittarius 
dSph was also removed from our sample because that is 
a clear case where tidal stripping has altered the lumi- 
nous component as well as the dark matter halo. This 
particular system will be studied in detail in subsequent 
papers of this series, which address the effects of tides on 
the results presented here. 

2.2. Modeling 

Our dwarf galaxies models assume the presence of two 
components in dynamical equilibrium: (i) a stellar com- 
ponent approximated by a King (1966) model, and (ii) a 
dark matter halo, which we approximate using an NFW 
profile (Navarro, Frenk & White 1996, 1997). We assume 
that the dark matter dominates the dynamics of the sys- 
tem, and that stars may be regarded as massless tracers 
of the potential. 

2.2.1. Luminous component 

The density profile of a King model may be written as 
(King 1962) 



f \ K 



cos 1 {x) 



where 
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(2) 



l + (r/r t ) 2 

Here p*(r > r t ) — 0; r k is an inner radial scale; r t is the 
King "tidal" radius; and K is an arbitrary normalizing 
constant. 

In the absence of rotation, the stellar kinematics is 
determined by the total gravitational potential, <&(r), 
theough Jeans' equations. In particular, the radial ve- 
locity dispersion of stars, ay, is given by 

r^p^ — dr, (3) 
dr 

(Binncy & Trcmainc 1987), where j3 is the velocity 
anisotropy. We shall assume hereafter that the stellar 
velocity distribution is isotropic (/3 = 0), which implies 
that the (observable) stellar line-of-sight velocity disper- 
sion, a p {R) is given by 

A ' Z{R)J R Vr 2 - R 2 ' 

where R is the projected radius and is the projected 
stellar density, 

m) = 2 P , = dr. (5) 



Vr 2 - R 2 
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We follow traditional convention and define the (pro- 
jected) core radius, R c , by the condition S(J? C ) = 
£(0)/2. The core radius, defined this way, depends on 
both of the King model parameters, and r t . For ex- 
ample, R c ~ 0.73 rk for r t /rk — 10. 

2.2.2. Dark matter component 

We assume that dark matter halos may be approxi- 
mated by NFW profiles (Navarro et al 1996, 1997). The 
density profile may be written as 



Pnfw 



M vir (r/r^jl + r/r,)- 2 
47rrf [ln(l + c) - c/(l + c)] ' 



(6) 



where M v ; r is the mass within the virial radius, r V i r , r s is 
the scale radius and c is the concentration (c = r v i T /r s ). 
The virial radius is defined so that the mean overdensity 
relative to the critical density is A v j r , 



(4/3)^ " Avir Pcrit ~ Avir 8irG 



(7) 



We follow Bryan & Norman (1998) and define the over- 
density by 



A vir (z) = 18tt 2 + 82 f(z) - 39 f{zf 



where 
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We assume throughout the paper a ACDM Universe 
by fixing the cosmological parameters to CIq = 0.3, 
£1a = 0.7, h = 0.7, n=l, and erg = 0.9, consistent with 
constraints from CMB measurements and galaxy cluster- 
ing (see Spergel et al 2006 and references therein). 

Finally, we note that an NFW profile is fully deter- 
mined by two characteristic parameters, such as the virial 
mass and the concentration. However, because the virial 
definitions adopted above depend on redshift, it is some- 
times preferable to characterize an NFW halo by the lo- 
cation of the circular velocity peak, (r max , V max ), where 
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(10) 

and r max ~ 2 r s . 

Although the two characterizations are equivalent, we 
shall adopt the latter in this paper, since it is independent 
of redshift and its parameters are more easily compared 
with observation. 

3. KING MODELS EMBEDDED IN NFW HALOS 
3.1. Spatial segregation and velocity dispersion 

Stars that follow a King model and inhabit an NFW 
halo have their kinematics dictated principally by the 
spatial segregation of stars relative to dark matter, which 
may be quantified by the ratio of King "core" radius to 
NFW scale radius, R c /r s . The more concentrated the 
stars are relative to the dark matter, the smaller the 
velocity of the stars will be relative to the characteristic 
circular velocity of the dark matter dominated potential. 

Fig. [T]shows the projected density profile of King mod- 
els as a function of R c /r s , and compares them to an NFW 
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Fig. 1. — Stellar surface density profile of King models with 
the same "concentration", Rt/R c = 10, and various degrees of 
spatial segregation relative to the dark matter. The segregation is 
measured by the ratio R c /r s , where r s is the scale radius of the 
NFW profile (shown by a dotted line in the figure) and R c is the 
stellar core radius, denoted with a dot in each profile. All radii 
have been scaled to the scale radius of the NFW profile. Units of 
surface density are arbitrary, since stars are assumed to contribute 
negligibly to the potential of the system. 




log( R/r s ) 

Fig. 2. — Stellar projected velocity dispersion profiles normalized 
to the halo maximum circular velocity for different values of R c /r s . 
The dotted line represents the projected velocity dispersion profile 
of the NFW halo. The solid circles indicate the core radius of the 
King model. Note that both the central velocity dispersion of the 
stars, as well as the shape of the velocity dispersion profile, depend 
on the degree of spatial segregation between stars and dark matter 
(quantified here by the ratio R c /r s ). 

profile. The radial units are scaled to the NFW scale ra- 
dius, r s , and the vertical units are arbitrary in this plot. 
The corresponding velocity dispersion profiles are shown 
in Fig. [2l Velocities have been scaled to the peak circular 
velocity of the NFW halo, Vmax, and show clearly the an- 
ticipated behaviour: the deeper the stars are embedded 
within the halo, the smaller the stellar velocities. For 
example, if R c ~ 0.1r s , the stellar central velocity dis- 
persion, (7p(0), is only about 40% of Vmax, but this value 



rises to 80% for R c ~ r s . There is also a weak depen- 
dence on the tidal radius, but this is minor, as discussed 
below. 

This implies that a "family" of NFW halos is consis- 
tent with a King model of given R c and er p (0) (sec Strigari 
et al 2006 for a similar argument applied to the Fornax 
dSph). The more embedded we assume the King model 
to be within the NFW halo, the more massive (i.e., higher 
Vmax) the halo must be in order to explain a given er p (0). 
This "King-NFW degeneracy" is illustrated in Figure [Sj 
where the thick lines show the correspondence between 
V max /a p (0) and r max /i? c . The three thick lines (almost 
indistinguishable from one another in this panel) corre- 
spond to three different values of the tidal-to-core radius 
ratio chosen for the King model, and confirm the result 
anticipated above regarding the weak dependence on R t 
of these results. 

Any NFW halo whose circular velocity peaks some- 
where along the thick curve in Figure [3] is, therefore, con- 
sistent with a King model of given cr p (0) and R c . A few 
examples of NFW halos belonging to this "family" are 
shown by the thin lines in Figure [3] Note that essentially 
all the circular velocity profiles of these halos, despite 
having very different masses and concentrations, cross 
each other at r w R c . This implies that the total enclosed 
mass within the core radius of a King model is robustly 
determined given our assumptions: in particular, we 
find V C (R C ) ~ 1.2ct p (0) and M{R C ) ~ 1.44 R c a p {0) 2 /G. 
Physically, this means that, although the total mass and 
radial extent of the halo are not well pinned down, the 
mass within the luminous region sampled by the tracers 
is. (See also Strigari et al. 2007 for a similar approach 
and conclusion.) 

A further constraint may be gleaned from Fig. [2j which 
shows that the shape of the stellar velocity dispersion pro- 
file also depends on the degree of segregation between 
stars and dark matter. For R c ^0Ar s (= 0.05 r max ) the 
velocity profile remains approximately flat well outside 
the core radius, and declines abruptly only at the "tidal" 
radius. On the other hand, less segregated King models 
show a steep velocity decline noticeable near the core ra- 
dius: for example, for R c ~ r s the velocity dispersion de- 
clines by roughly 70% at the core radius from the central 
value. The velocity dispersion profiles of all dwarfs for 
which such data are available show little sign of declin- 
ing outside R c (see references in § 12. 1[) . In the context of 
our modeling, this suggests that the stellar component 
is deeply embedded within its parent CDM halo (i.e., 
Rc ~ ^max), an issue to which we return below. 

3.2. Application to Milky Way dwarfs 

One way of breaking the degeneracy illustrated in Fig- 
ure [3] is to appeal to the results of cosmological N-body 
simulations. These show that there is a strong correla- 
tion between the mass and concentration of a CDM halo 
or, equivalently, between r max and V max . As discussed 
by NFW, this correlation arises because the characteris- 
tic density of a halo is proportional to the density of the 
Universe at the time of its assembly. In practice, and 
for galaxy-sized systems, the mass-density dependence is 
quite weak, implying that r max is roughly proportional 
to Vmax. This correlation is now well established, and a 
number of authors provide simple formulae to compute 
it once the cosmological parameters are specified (NFW, 
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Fig. 3. — The King-NFW degeneracy. The thick lines show 
the halo peak circular velocity, I4nax, in units of the central ve- 
locity dispersion, and the radius of the peak, r max , in units of 
the King model core radius. Solid, dotted and dashed lines de- 
note different King-model concentrations (Rt/ Rc)- Any NFW halo 
whose circular velocity peaks along this curve is consistent with 
the King model structure and kinematics. A few of these NFW 
models are shown for illustration by the thin curves. Note that 
all these NFW models cross each other at approximately R ~ R c 
and V c — 1.2(j p (0). This implies that, given our assumptions, the 
mass within the core radius of the stellar component is robustly 
constrained to be M(R C ) ~ 1.44 R c ct|(0)/G. 

Eke, Navarro & Steinmetz 2001, Bullock et al 2001). 

The relation between V^ nax and r max implies that, of all 
NFW halos in the "family" of models allowed for a given 
dwarf by the degeneracy illustrated in Figure El a single 
one will be consistent with the parameters expected in a 
given cosmogony. This is shown in Figure [H where the 
curved lines show, in different panels and for each of the 
eight MW dSphs, the King-NFW "degeneracy" relation. 
The straight lines in each panel delineate the V m&x -r max 
relation consistent with the ACDM cosmogony. The set 
of three straight lines correspond to NFW halos identified 
at z = 0, 1, and 2, respectively. We shall hereafter ignore 
the relatively small difference between these curves and 
refer, for simplicity, all of our results to z = 0. 

In each panel a solid dot indicates, for reference, the 
core radius and the central velocity dispersion of each 
galaxy, as listed in Table [T] All of these points lie well 
to the left of the cosmological relations, confirming our 
earlier suggestion that the luminous components are sig- 
nificantly segregated within their dark halos (i.e., they 
are substantially smaller for given characteristic veloc- 

ity). 

The intersection between the King-NFW degeneracy 
and the cosmological relation is marked by red symbols 
for each dwarf. This indicates the parameters of the 
NFW halo model that is consistent with the ACDM cos- 
mogony and, at the same time, matches the kinematics 
and structure of each dwarf galaxy. The circular velocity 
profiles of the eight NFW halos satisfying these criteria 
are shown in Fig. [5l and labelled from top to bottom in 
order of decreasing halo mass. 

This ranking shows interesting peculiarities. For ex- 
ample, it shows that the peak circular velocities of dwarf 
halos vary from ~ 17 to ~ 35 km/s, corresponding to 
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Fig. 4. — Each panel shows, for the 8 Milky Way dwarfs in 
our sample, the King-NFW degeneracy (curved line), as well as 
the predictions for ACDM cosmogony (set of straight lines). Only 
NFW halos at the intersection of both sets of curves (marked by 
a red symbol) are consistent with both cosmological constraints 
and the structure and kinematics of the dwarfs. The three sets of 
cosmological curves correspond to NFW halos identified at vari- 
ous redshifts; we adopt the 2 = models here, but note that our 
conclusions are unlikely to be severely affected by this choice. 

a spread of about 8 in mass, much narrower than the 
factor of ~ 70 spanned by dSph luminosities. Note as 
well that, as anticipated in §1, halo mass is not mono- 
tonically related to luminosity. Intriguingly, Draco, one 
of the faintest dwarfs in our sample, is assigned the most 
massive halo, whereas Fornax, despite being 70 times 
brighter, is assigned a halo 5 times less massive. 

The lack of correlation between luminosity and halo 
mass is shown explicitly in Figure [51 where we plot, as a 
function of total visual magnitude, the virial mass of the 
halo (red symbols at the top) , as well as the mass within 
the core radius of each dwarf (blue symbols at bottom). 
Symbols are the same as those used in Fig. [5] to denote 
different galaxies. 

Note that the total mass within the core radius — 
a fairly robust measure according to our discussion of 
Fig. [3] is approximately independent of luminosity 4 . 
This implies that the average density of dark matter will 
be higher in physically smaller systems such as Draco 
than in more extended dwarfs such as Fornax. Why does 
this matter? Because more massive halos are denser than 
less massive ones at all radii in the CDM cosmogony. In- 
deed, note that the circular velocity profiles of the NFW 
halos shown in Figure [5] do not cross, which means that 
measuring the halo circular velocity (or mass) at any ra- 
dius leads to a well-defined estimate of the total mass of 
the halo. Because Draco, despite being faint, has a cir- 
cular velocity comparable to Fornax's at a much smaller 
radius, it requires a denser, and therefore more massive, 
halo to satisfy the observational constraints. 

Of the eight dwarfs, the lowest halo mass corresponds 
to Sextans, whose relatively large radius and small veloc- 
ity dispersion is inconsistent with a very massive CDM 

4 This is not a new result, and is consistent with Mateo's (1998) 
conclusion that simple dynamical mass estimates of dwarf galaxies 
are independent of luminosity. 
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halo. 

3.3. Effects of tidal stripping 

The results obtained in the previous section assume 
that the structure of dark matter halos is well approx- 
imated by an NFW profile. Although this assumption 
may be appropriate for isolated halos, it is unlikely to 
hold in detail for "substructure" halos of dSphs orbit- 
ing within the main halo of the Milky Way. Recently, 
Stoehr et al (2002), Hayashi et al (2003) and Kazantzidis 
et al (2004) have examined the modifications undergone 
by an NFW halo as it is tidally stripped inside a more 
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massive system. Stripping affects principally the outer 
regions of the halo, and therefore deeply embedded stel- 
lar structures may survive unscathed the removal of large 
fractions of their halos. 

This is illustrated in Figure El which shows the circu- 
lar velocity profile of an NFW halo, after being tidally 
stripped of 50%, 75%, and 90% of its mass, respectively. 
The profiles are taken from N-body simulations of NFW 
halos orbiting within the potential of a much larger sys- 
tem, as in Hayashi et al (2003). Figure [7] indicates 
that NFW halos must lose at least 90% of their original 
mass before regions within r^50.1r max are significantly 
affected. Interestingly, this is precisely the region that 
our analysis suggests is populated by the stars of MW 
dwarfs. According to Fig. [5] the average i? c /^max for all 
MW dwarfs is 0.054 and they all have, with the possible 
exception of Sextans, R c ^ 0.1 r max . 

Tidal stripping is clearly affecting the structure of the 
Sagittarius dSph (which is thus not in the sample consid- 
ered here), but the evidence is less clear-cut in the case 
of other dwarfs, despite a rich literature on the topic. 
Much of this work is based on detecting stars associated 
with a dwarf beyond the King "tidal radius" , or on the 
interpretation of "breaks" in the surface density profile in 
the outer regions of a dwarf (e.g. Ursa Minor, Martmez- 
Delgado et al. 2001, Carina, Majewski et al. 2005). 

The presence of a break in the outer profile, how- 
ever, does not imply that the system is necessarily losing 
stars. Such feature might just be intrinsic to the dSph, 
as shown, for example, by the multiple-component mod- 
els presented in McConnachie et al. (2006). These au- 
thors show that systems with multiple stellar components 
(such as Andromeda II) may have overall density profiles 
that resemble those of dwarf galaxies with outer "tidal" 
breaks. The presence of well-defined and kinematically- 
distinct populations of stars within some dSphs (Sculp- 
tor, Tolstoy et al. 2004; Fornax, Battaglia et al. 2006; 
Canis Venatici, Ibata et al. 2006), indeed, argues against 
tides as a major driver in the evolution of a dSph, since 
tidal stirring would tend to uniformize such distinctions. 

On similar grounds, one cannot rule out that other fea- 
tures in the light profile usually ascribed to tides, such 
as the presence of lumps (UMi, Olszewski & Aaronson 
1985), shells (Fornax, Coleman et al. 2005) and asphcri- 
cal isopleths (UMi, Martmez-Delgado 2001) may actually 
reflect complexities in the formation process and/or the 
internal evolution of dwarf galaxies. 

The lack of overwhelming evidence for ongoing strip- 
ping of stars from most dSphs may thus be taken to imply 
that dwarf halos have retained at least 10% of their orig- 
inal mass at present. This is important, because even for 
such sizeable loss, the peak circular velocity of the halo 
is barely affected. 

The thin line in Figure [7] tracks the position of the 
circular velocity peak of an NFW halo as it is tidally 
stripped, and shows that halos that have lost 75% of their 
mass to stripping see only a 10% decrease in Kn ax . Even 
after losing 90% of its mass, a halo sees its U max reduced 
by only ~ 30% 5 . We conclude, therefore, that, unless 
dSphs have suffered catastrophic tidal losses, our esti- 

5 The location of the peak, r max , on the other hand, is more 
significantly affected, and shifts inwards by almost a factor of ~ 3 
after a halo loses 90% of its mass. 
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Fig. 7. — NFW halo circular velocity, normalized to the location 
of the peak. The core radius of each of the MW dwarf galaxies is 
shown. Note that, except for Sextans, in all cases the core radius 
lies within 0.1r max . Dotted lines show the changes imposed by 
tidal stripping, as the original NFW halo loses 50%, 75%, and 90% 
of its mass. Note that most dwarfs are sufficiently embedded within 
their halos that they would only be affected if the halo lost more 
than ~ 90% of their original mass. 

mates of the peak circular velocities for the MW dSphs 
(Fig. [5]) should be relatively robust, even if the amount 
of dark matter beyond the luminous radius remains un- 
known. 

We emphasize, however, that tides do pose a num- 
ber of interesting questions that our modeling fails to 
address. For example, do King profiles survive strong 
tidal stripping?. Do we expect the presence of "extra- 
tidal" stars beyond a break radius if the luminous com- 
ponent is still surrounded by dark matter? These issues 
are best addressed by self-consistent N-body simulation 
of the tidally-driven evolution of multi-component dSph 
models. This is beyond the scope of the present paper, 
but will be dealt with in a forthcoming paper of our se- 
ries. 

3.4. Application to the "satellite crisis" 

The analysis of the preceding subsections suggests that 
the 8 MW dwarfs chosen for our analysis inhabit halos 
with peak circular velocities in the range 17-35 km/s, 
a factor of ~ 3 higher than their central (stellar) ve- 
locity dispersions. As discussed by Stoehr et al (2002), 
Hayashi et al (2003) and Kazantzidis et al (2004), this 
correction is enough to alleviate substantially the "satel- 
lite crisis" highlighted by Klypin et al (1999) and Moore 
et al (1999). This is illustrated in Fig. [51 which com- 
pares the estimates of the peak circular velocities ob- 
tained in this contribution (full circles) with those that 
adopt mass-follows-light models (full triangles). 

In rough terms, N-body simulations show that there 
are typically about 20-30 substructure halos with peak 
circular velocities exceeding 10% of the main halo's virial 
velocities. That is shown in Fig.[5]from the results of Dic- 
mand, Kuhlen & Madau (2006) for two different Milky 
Way- like halos (dotted and dahsed lines). According 
to our analysis, there are about a dozen dwarfs with 
U m ax ~ 17 km/s (adding to our sample brighter satel- 
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Fig. 8. — Cumulative number of substructures as a function of 
maximum circular velocity for two Milky Way-sized ACDM halos 
(solid and dashed lines, data from Diemand et al. 2006). Trian- 
gles approximate the maximum circular velocity of the Milky Way 
satellites as V ma x — \/2a p (ff) (Klypin et al. 1999, Moore et al. 
1999) whereas circles show the estimates of V max obtained in this 
contribution. For the LMC and SMC we assume that V max equals 
the maximum rotation velocity of the HI disk (Kim et al 1998, 
Stanimirovic et al 2004). Open symbols denote systems with mea- 
sured and lacking velocity dispersion measurements, respectively. 
Note that our estimates considerably alleviate the "satellite crisis" 
for satellites with peak circular velocities exceeding 10-20 km/s. 
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Fig. 9. — Comparison of the surface density profile of pairs of 
dwarf galaxies of similar luminosity, one orbiting around M31 and 
the other around the Milky Way. On average M31 dwarfs are larger 
and have lower surface brightness than their MW counterparts. 

lites such as the Magellanic Clouds) . This velocity corre- 
sponds to between 8 and 11% of the MW virial velocity, 
assuming, as seems likely, that the latter is in the range 
150-220 km/s. This factor of ~ 2-3 discrepancy between 
the number of massive dark matter substructures and lu- 
minous satellites does not appear extravagant given the 
uncertainties. 

Furthermore, the discrepancy may disappear alto- 
gether if (as the scenario outlined in §[1] would suggest) 
the ultra-faint MW dwarfs identified in SDSS data (open 
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FlG. 10. — Upper panel: Core radii of Milky Way (filled circles) 
and M31 (filled squares) dwarf galaxies as a function of their ab- 
solute magnitude. The dashed lines are intended to guide the eye 
and to highlight the difference in size between the M31 and MW 
dwarf populations. Bottom panel: Measured (filled symbols) and 
predicted (open symbols) central velocity dispersions if the Milky 
Way and M31 dSph's are embedded in similar dark matter halos. 
See text for details. 

symbols in Fig. [5]) turn out to inhabit halos of masses 
comparable to those of the 8 dSphs we consider here. 
Since these newly discovered dwarfs have, on average, 
similar physical size to the dSphs in Table 1, we expect 
that their velocity dispersions will also be comparable 
and in the range ~ 6-10 km/s, a prediction of this CDM- 
motivated modeling that should be testable in the near 
future. 

Indeed, further progress on this subject seems immi- 
nent, as observational campaigns secure velocity disper- 
sions for many of the newly discovered dwarfs and ex- 
tend available data to allow for velocity dispersion pro- 
files to be measured for more systems (see, e.g., the re- 
cent preprint by Simon & Geha 2007). One also expects 
that numerical simulations will improve to the point that 
a reliable direct estimate of the number of substructure 
halos satisfying the M(R C ) constraints shown in Fig. [6] 
will be possible, obviating the need to extrapolate results 
out to the peak of the halo circular velocity curve. The 
latter goal requires simulations able to resolve convinc- 
ingly the inner ~ 100 pc of substructure halos, an order 
of magnitude improvement over simulations published so 
far. This is a numerical challenge that will likely be met 
soon but that will nevertheless require the investment of 
massive computational resources (see the recent preprint 
by Kuhlen et al 2007 for a report of recent progress on 
this issue.) 

4. MILKY WAY VS. M31 DWARFS 

Photometric studies of the dwarf galaxies in the Local 
Group (Irwin & Hatzidimitriou 1995, McConnachie & 
Irwin 2006) have shown that, at fixed luminosity, M31 
dwarfs are considerably more extended than their Milky 
Way counterparts (see Table [I}. We illustrate this in 
Fig. where King model fits to the surface brightness 
profiles of three pairs of dwarfs of similar luminosity are 
shown. Each pair consists of one M31 satellite and one 
MW satellite, respectively. Clearly, the M31 dwarfs are 
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systematically different: they are about 1.5 mag fainter 
in central surface brightness and about a factor of two 
larger in radial extent compared with MW dwarfs (see 
top panel of Fig. [TUj) . 

The origin of this difference is unclear, but it offers 
nonetheless a way of assessing the general validity of the 
modeling proposed in the previous sections. Indeed, un- 
der the plausible assumption that the dark matter halos 
of dwarfs in M31 and in the Milky Way are not sys- 
tematically dissimilar, the arguments of § 13. II imply that 
the difference in size of M31 dwarfs should be reflected 
in their kinematics through significantly higher velocity 
dispersions. Note that this prediction depends explicitly 
on the presence of dark matter halos whose mass and ex- 
tent are largely independent from the stars', as implicitly 
assumed in our modeling so far. Other (also plausible?) 
assumptions, such as mass-traces-light, would lead to the 
opposite trend, and would predict systematically lower 
velocity dispersions. 

In order to estimate quantitatively the velocity dis- 
persion of M31 dwarfs we need to know how spatially 
segregated the stars are from the dark matter in these 
systems. Once this is fixed by, for example, assuming a 
value for the ratio R c /r max , we can compute r max for each 
dwarf and, from that, infer Vrn ax assuming the cosmolog- 
ical relation between these parameters. Since the ratio 
i? c /r max also fixes V max /a p (0) according to the King- 
NFW degeneracy shown in Fig. [3] an estimate can then 
be made of the central velocity dispersion, o~ p (0), for each 
M31 dwarf. 

We have followed this procedure, assuming i? c /r max = 
0.05, consistent with the average segregation measure de- 
rived for MW dwarfs (see Fig. [7]) , as well as with the flat 
o-p(R) constraint mentioned in § 13.11 This leads to the 
estimates of the velocity dispersions of M31 dwarfs plot- 
ted in the bottom panel of Fig. [10] (open symbols). Error 
bars show the variation induced in the estimate by allow- 
ing i? c /r max to vary between 0.025 and 0.1. Clearly, the 



velocity dispersion of M31 dSphs is significantly larger 
(by about a factor of ~ 2) than that of their MW coun- 
terparts, as shown by the shift between solid circles and 
open symbols in Fig. 1101 

We emphasize that these "predictions" for the velocity 
dispersions of M31 dwarfs should be taken with caution: 
indeed, the open symbols in Fig. [TUJ are best interpreted 
as indicative of the relative shift in the velocity dispersion 
of the population of M31 dwarfs rather than as separate 
individual predictions. 

So far, the few available observations are consistent 
with the predicted trend. Of the three dwarfs outside 
the MW for which velocity dispersions have been pub- 
lished (see, Table[JJ), the most reliable is Cetus and, at 17 
km/s, has the highest o~ p in the sample, consistent with 
the predicted trend (our naive prediction yields H-4T^ 
km/s). Likewise, the higher of the two estimates derived 
by Chapman et al (2005, 12 km/s) for And IX is also 
consistent with this trend, although the very few stars 
in this galaxy (as well as in And II) with measured ve- 
locities suggest that it would be premature to draw firm 
conclusions until the data improves substantially. 

It is clear, however, that should future observations 
fail to confirm the predicted trend, one or many of our 
assumptions would need to be revised. As an illustra- 
tive example, we consider how our interpretation would 
change if the velocity dispersion of M31 dwarfs were to 
show no systematic departure from that of their MW 
counterparts (i.e., if the open symbols were found to align 
with the dashed line in the bottom panel of Fig. [TO]) . 
In this case, one would be forced to conclude that M31 
and MW dwarfs of similar luminosity inhabit systemati- 
cally different halos. In particular, M31 dwarf halo hosts 
would need to be significantly less massive (i.e., lower 
Knax and smaller r max ) in order to accommodate their 
larger physical size but similar kinematics. 

This is shown in Fig. [TTJ and imply that essentially all 
M31 dwarfs considered here would have V max ^> 15 km/s 
(with the possible exception of And VII). Since, in all 
likelihood, M31 inhabits a halo at least as massive as the 
Milky Way's, it should have at least as many massive 
substructures as our own Galaxy, and one would need 
to explain why those substructures have not "lit up" in 
M31 as they have in the Milky Way. We consider this a 
compelling argument against this interpretation and in 
support of higher velocity dispersions for M31 dwarfs as 
the most natural prediction of CDM-motivated models 
of the Local Group dSphs. 

One last possibility should be mentioned; namely, that 
the structural differences between M31 and MW dwarfs 
reflect differences in the dynamical evolution driven by 
tides in M31 and the Galaxy. Although we have ar- 
gued in §3.3 that this is unlikely, a definitive assessment 
requires high-resolution, self-consistent N-body simula- 
tions of how tides affect multiple-component models of 
dSph galaxies. We plan to address this issue in the foll- 
woing paper of this series (Peharrubia, Navarro & Mc- 
Connachie, in preparation). 

5. SUMMARY 

We have considered in this paper the observational 
constraints placed on the mass and spatial extent of 
cold dark matter halos surrounding Local Group dwarf 
spheroidal galaxies. Assuming that the luminous com- 
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ponent may be approximated by King models and that 
dark halos follow the NFW mass profile, we conclude that 
estimates of the halo mass of a dSph depend principally 
on the degree of spatial segregation between stars and 
dark matter: the more embedded the stars within the 
dark halo the more massive the halo must be in order to 
explain the observed stellar kinematics. This degeneracy 
may be broken by appealing to the results of cosmological 
N-body simulations — which indicate a strong correlation 
between mass and size of cold dark matter halos. The 
procedure results in reasonably robust estimates of the 
dark halo mass of a dSph and of its peak circular velocity. 

Our analysis indicates that the mass within the lumi- 
nous radius of the dwarf is well constrained by the veloc- 
ity dispersion of the stars, and that it is approximately 
independent of the luminosity of the system. Within 
the context of cosmologically motivated CDM halos, this 
implies that dSphs are surrounded by halos with peak 
circular velocities a factor of ~ 3 times larger than their 
stellar velocity dispersion. These results are consistent 
with previous work, and substantially alleviate the "miss- 
ing satellites" problem, particularly if most of the newly 
discovered ultra-faint satellites are confirmed to be bona- 
fide dwarf galaxies akin to the more luminous, classical 
systems analyzed here. 

We also find that stars in dSphs are deeply embedded 
within their dark matter halos, with core radii typically 
of order 5% the radius where the original halo circular 
velocity peaks. This yields velocity dispersion profiles 
that are nearly flat out to the nominal King tidal radius, 
as observed for most of the Milky Way dSphs with suit- 
able kinematical data. The deep segregation of the sellar 
component also implies that the stars of dwarf galaxies 
are fairly resilient to tidal disruption; a halo would need 



to lose more than 90% of its mass before stars begin to 
be affected. Even in this case, the peak circular velocity 
of a halo would only drop by ~ 30%, suggesting that our 
U max estimates are robust even if halos have been heav- 
ily (but perhaps not catastrophically) affected by tidal 
stripping. 

Applied to the dSph population of M31, this analysis 
suggests that the systematic difference in size between 
M31 dwarfs and their MW counterparts should be re- 
flected in their kinematics. Under the plausible assump- 
tion that substructure halos in M31 are similar to those 
of the Milky Way, we conclude that, as a population, the 
velocity dispersion of M31 dwarfs should be ~ 50-100% 
higher than that of dwarfs in our own Galaxy. 

This prediction is likely to be validated (or challenged) 
soon by the results of the many ongoing observational 
projects devoted to obtaining accurate spectra and ra- 
dial velocities for stars in these galaxies. Whether Local 
Group dwarfs conform with the expectations of the pre- 
vailing CDM paradigm or throw it a further gauntlet we 
should know in the near future. 
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